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Spin foams are candidate state-sum models for transition amplitudes in quantum gravity. An active 
research subject is to identify the possible divergences of spin foam models, or alternatively to show 
that models are finite. We will discuss in detail the (non-occurrence of) divergences in the Barrett- 
Crane model, formulated as an integral of delta function weights only. We will furthermore present 
a simple method to estimate the divergence degree of the so-called bubbles for general spin foam 
models. 

Divergences in spin foams are expected to be related to the existence of gauge symmetries (dif- 
. . . feomorphisms) . Thus we have to conclude that such gauge symmetries are not (fully) present in the 

CO I model we consider. But we will identify a class of gauge symmetries which occur at special solutions 

of the equations imposed by the delta function weights. This situation is surprisingly similar to 
^^ ' the case of broken diffeomorphism symmetries in discrete gravity, which are present around flat 

Cn ' solutions. We introduce a method to derive (Ward-identity-like) equations for the vertex amplitude 

}_( , of the model in the case of broken gauge symmetries. 
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I. INTRODUCTION 
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I ' Spin foam models arose as a path integral approach to quantum gravity yj. One of the first specific 4D gravity 

wp.! models constructed was the Barrett-Crane (BC) model 3, a more modern alteration of which is 0]. The EPRL/FK 
models Q were also recently proposed due to arguments that the Barrett-Crane model fails certain tests to be a 
viable model for gravity [5|. 

All models underly a unifying construction principle [6[, which in our view is quite accessible in the holonomy 
representation [7|. There the common starting point is the use of distributions (delta functions) as weights for the 
1^^ . faces of the 2-complex on which spin foam models are defined. These delta function face weights are altered by edge 
\^ ■ functions convoluted into the face weights. These edge weights are in general also distributional. For the BC model 
\^ , these are again delta functions (for a specific choice of what is termed edge weight or measure factors), and this will 

allow us to evaluate quite explicitly certain configurations, 
f-k I Given that the models involve distributional objects one has to worry about divergences and the question arises 

^r^ . whether there is need for a regularization. This is an actively studied issue jsl-liot relevant for the definition and 
behaviour of group field theories |ll| . for which spin foam models provide the Feynman amplitudes, and the regular- 
ization of spin foams in itself. 

In this work we will consider in detail the possible divergences that can occur in the Barrett-Crane model (with a 

^^ , specific choice of edge and face weights, which has not been considered before). We will also present a simple method 

H ■ to estimate the occurrences of (single bubble) divergences in general spin foam models (including EPRL/FK). This will 

- - -' show explicitly how the choice of face and edge weight factors influences the divergence properties of the models. Here 

it will turn out that one has to choose between the invariance of the model under certain edge and face subdivisions 

and convergence. If one wants to avoid single bubble divergences (divergences related to diffeomorphism symmetry 

would be expected from multiple bubble configurations) one might also take convergence considerations into account 

in the determination of these factors [l2| - KL4J . 

Apart from the possible need for regularization there is another strong motivation to study the divergence structure 
of spin foams. This is the relation between divergences, the redundancies of delta-distributions and diffeomorphism 
symmetry Il5tjl9ll. and the proposal that spin foams act as projectors onto the Hamiltonian and diffeomorphism 
constraints [2Cl |2l| . As diffeomorphism symmetry leads to non-compact orbits (for vanishing cosmological constants) 
one would expect that an anomaly free implementation [Tj] of these symmetries would lead to divergences due to the 
integrations of amplitudes over the non-compact orbits(on which these amplitudes are constant). Indeed this relation 
is well understood for the 3D Ponzano-Regge model [ISl . Il7| , where the diffeomorphism symmetry is implemented as 
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a translation symmetry on the vertices of the triangulation. This symmetry also allows the derivation of recursion 
relations [22[ that can be related to the Hamiltonian [23, [2J] . 

The issue is more involved for the 4D (gravitational) models. Spin foams can be seen as a discretization of the 
path integral. Thus the question arises whether diffeomorphism symmetry is preserved under discretization even on 
the classical level. The single simplex amplitudes of the spin foam models approach the Regge action in the large 
j limit [2^ ■ The 4D Regge action, as a discretization of the Einstein Hilbert action does however break in general 
diffeomorphism symmetry (as opposed to the 3D Regge action) [26[. This statement has however its exceptions 
[19l |26|: The subdivision of a 4-simplex into five simplices by placing an inner vertex in the inside of the initial 
vertex leads to a configuration with vertex translation symmetry. This also leads to (classical) first class Hamiltonian 
constraints J26l - l28l | for four-valent vertices in a triangulated three-dimensional hypersurface. Thus one could expect 
divergences for the spin foam models at least for these configurations. The corresponding symmetry (for instance in 
the form of redundant delta-functions) could then be used to derive recursion relations and a quantum Hamiltonian. 

We will however argue that divergences for the BC model only appear for very special configurations, that first of 
all have to include two-valent faces (i. e. faces with only two edges) and furthermore have to combine these two-valent 
faces in a specific way. (Indeed a gauge symmetry can be found for the case that two two-valent faces are glued onto 
each other. This is the only configuration for which we found a divergence.) Thus we do not expect a full gauge 
symmetry related to the subdivision of a simplex, which does not involve two-valent faces. 

The subdivision of a simplex corresponds to a situation where aljj (classical) solutions are connected by a gauge 
symmetry. Another case is the occurrence of special (i.e. fiat) solutions on more general triangulations. The Hessian 
around these solutions will feature null modes which signifies the existence of gauge symmetries around these special 
solutions. 

Such symmetries have not been discussed for spin foams so far. Here we will consider an analogue situation for spin 
foams, that is analyze (gauge) symmetries that occur around special solutions. This is also the reason for considering 
mostly the BC model in this paper, as (with our choice of edge and face weights) it can be rewritten as the integral 
over a space of flat connections, i.e. a partition function with only delta function weights. The special solutions are 
special points in the space of fiat connections. In this work we will present a method to derive recursion relations for 
the vertex amplitude of the BC model, the lOj symbol, which are derived from these special solutions. 

The structure of the paper is as follows. After introducing the Barrett-Crane model with our choice of face and edge 
weight factors in section |TT1 we will switch to a group integral formulation in section IIIIl This will introduce effective 
face weights which capture the possible divergences of spin foam models. For the BC model these effective face weights 
can be evaluated explicitly and we find them finite for faces with more than two edges (modulo divergences which 
occur on measure zero sets). After considering the square of such effective weights, showing that the measure zero 
set divergences do not matter in this case, we continue with a discussion of so-called bubble divergences in sections 
IIVI The methods used there can be generalized to other models as well. As we will show in section HV CI this allows 
for a simple estimate of possible divergences occurring in spin foam models. We then discuss the multiple bubble 
case, in particular the 4-dipole configuration, important for group field theories, in section |Vl To this end we will 
reformulate the partition function as an integral over a space of fiat connection. This technique will be essential for 
the consideration of gauge symmetries, which for the BC model occur around special solutions. We will use this 
technique in order to discuss these gauge symmetries in section IVIII and use these symmetries in order to derive 
recursion relations for the lOj symbols. We close with a discussion and outlook in section IVIIII The appendix \X\ 
includes some necessary basics on the group SU(2). 

II. THE BARRETT-CRANE MODEL 
A. Presentation 

In this section we will shortly introduce the model we will be considering in the rest of the paper. 

Let r be a two-complex. Ti its set of i-cells. We call 0-cells vertices, l-ccUs lines and 2-cclls faces. The spin 
representation of the Barrett-Crane (BC) model is as follows. A state is an assignment of spins {jf £ ]N/2}/gr2 to 
faces. To simplify, assume F is the 2-skeleton of the dual to a four-dimensional triangulation. Each vertex is dual to a 
4-simplex and each line to a tetrahedron. A 4-simplex has five boundary tetrahedra, so that a vertex in F has degree 



^ Here we assume that all classical solutions for this case are flat. There could be some special solutions which correspond to discretization 
artifacts however [2611. 



5. The faces of F are dual to triangles. Consider a vertex in T, and denote the ineoming lines a — 1, ... ,5. There are 
ten faces which are identified as the pairs of lines, (a, 6) for 1 < a < 6 < 5 (and corresponding to the ten triangles of 
the dual 4-simplex). Each vertex receives a weight, known as the lOj-symbol, labeled by the ten spins of the faces 
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{niab} = / n "^'^^ n 

^SU(2)= ,^1 l<a<&<5 



XjatiK^h) = 




(1) 



Here Xj is the SU(2) character in the representation of spin j, Xj( 



jei 



Em=-,e™'' = sin(dj6')/sin6', with n 



a normalized 3- vector, a ~ {o'x,'^y,'^z) the vector formed by the Pauli matrices and the notation dj — 2j + 1. The 
partition function of the model is a state-sum, 



^Bc = E n ^/ n ^^ n iio-?}- 



(2) 



veFo 



Several choices can be found in the literature for the measures on faces Af and lines Ag . We will stick to the following 
choice. 
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B. Finiteness 



The above choice of measure makes the model quite convergent. This can be found as follows. There are as many 
sums as faces in F. First we collect the factors djj for each sum. The face measure brings dj for any face. For a 
face with n vertices, dispatching the line measure A^ on the faces yields a factor rf~". Next we have to deal with the 
lOj-symbols. Previous studies strongly suggest a bound of the type 



\{lOj}\<K l[d 



3f 



(4) 



for some positive a. In [29J, numerical evidence gives a = 1/5. Therefore, 



i^Bci< E nc^ni{ioj}i^^'n 
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As a result, if all faces have at least 7^/ > 3 vertices, the partition function is finite. However, this result is far from 
satisfying as it does not provide any insight for quantum gravity. Here are some questions that it leaves unanswered. 

• The above arguments fails in the presence of faces with two vertices only. Such faces appear in generalized 
triangulations, in particular in the melonic sector which dominates group field theories. We would like to be 
able to treat them and understand why they are more likely to bring divergences from the quantum gravity 
point of view. 

• Any spin foam model can be made finite (or arbitrarily divergent) by adding negative (positive) powers of dj 
in the face and line measures. Several such different versions of the Barrett-Crane model exist [3, |30|. We will 
have to justify our choice ^ . 

• We would like to formulate the model on arbitrary two-complexes, not only those 2-skeleta dual to regular 
triangulations. This is necessary to understand bubble divergences (when only part of F is taken into account). 

• In the BF model, the Wilson loops are constrained to be trivial, and typical divergences come from redundancies 
in the set of constraints. When turning it to the BC model, are there constraints left? Is the above finiteness 
result related to absence of constraint redundancies? 



• These redundancies in the BF model are associated to the existence of gauge symmetries which survive on the 
lattice (and are not gauge-fixed). Does the above finiteness result imply the absence of such gauge symmetries 
in the BC model? 

We will answer these questions in the remaining of the paper, through a group integral formulation of the model, 
where sums over spins are traded for group integralqj. 

III. GROUP INTEGRAL FORMULATION 
A. Effective face weights 

There exists a group integral representation of this model, which comes out quite naturally as a derivation of the 



model from a path integral for discretized general relativity [31[. To each pair line- vertex (e, u) (where ti is a vertex 
of e) we associate an SU(2) element hg^ (they equivalcntly live on half-lines). They have the following geometric 
meaning. Each 4-simplex of the triangulation (dual to a vertex) is flat and equipped with a local Euclidean frame in 
]R^. Each boundary tetrahedron, which is dual to a line adjacent to the vertex, spans a 3-dimensional subspace of 
]R^, determined by the normal to it in the frame of the simplex, denoted 7Ve„ (with notations on F). This is a unit 
vector, hence an element of the 3-sphere. Using the isomorphism between the 3-sphere and SU(2), the normal N^v is 
representee^ as the element h^v of SU(2). 

The interest of such a representation of the normal is that geometric quantities can be expressed through the group 
law. The dihedral angle between two tetrahedra 61,62 in a 4-simplex v is the scalar product between their normals 
which can be written 

cos9ej^e2 = -tr h~^\he2v, (6) 

in the fundamental matrix representation of SU(2). 

A face / G F2 can have an arbitrary number n of lines and vertices on its boundary. We divide the face into wedges, 
where a wedge is a pair (/, w) or equivalcntly a pair of half-lines attached to v. There are n wedges Wif, . . . ,Wnf 
around /. The product of group elements entering ^ is canonically associated to a wedge w and to simplify notations 
we write 

The class angle of h^u is therefore the dihedral angle between the tetrahedra dual to the half-lines. We can now 
re- write the partition function as 

Zbc= / J]^ rf/iei, JJ a;(/iu,i^.,...,/ii„„^.), (8) 

{e,v) /era 

The function uj has n arguments and is called the effective face weight. It is fully determined by our choice of measures 
Af,Ae to be 

^ n 

L0(hi, . . . ,hn) = I TTc^7i S{hi J2h2j2^ ■■■-fnhnJn^), (9) 

where 6 is the Dirac delta over SU(2). 

To prove the equivalence between the spin representation ([2]) and the group representation ((8)) of the partition 
function, we start by expanding w onto the SU(2) modes, using S{h) = J2je-!M/2'^J^jW- ^ ^^^ more formula are 
needed in order to integrate the elements 7i in ([9]) . They are given in the Appendix [X] The character of products of 
group elements expands onto the Wigner matrices D'^^'> as follows Xi(^i^2) = 'J2m n=-i Dnin{hi)Dnm{h2)- Combining 
this with the orthogonality of the matrix elements of the Wigner matrices (|A4p . we get 

w(/ii,...,/i„) = 5]d2^"x,(/ii)x,(/i2)---X,(/i„), (10) 



^ Not surprisingly, integrals are easier to evaluate than sums. 

^ The map is N^^ = tr(/iei;0'^), for / = 0, 1, 2, 3, cr" = I, and a' the Pauli matrices. 



then 



^Bc = E / n ^'^^ n 4""' n xw(/^«.)- (n) 

{jf} {e,v) /6r2 wCf 

The factor d^ gives the face measure Af, the factors d- ^ are associated to pairs (e, /) and can be re-arranged as the 
hne measure A^. Finally the product of characters over faces can be re-organized as a product over vertices, yielding 
©• 

A first outcome of this formulation is that it directly makes sense on arbitrary two-complexes, not necessarily dual 
to a regular triangulation. 

The effective weight uj is obviously well-defined as a distribution. Let us integrate uj with some regular test function, 

/n ^ n 

2 — 1 Z— 1 

/n— 1 n 

J]^ dh^ Jl d7, -0(^1, • ■ • , K^l,ln^ {ln^lh~-lln-l) ' ' ' {nK^lI^) In) (12) 

4=1 i=l 

/n— 1 n — 1 

n ^5» n '^^i '^{ll^am, ■■■, 7r7-l5n-l7n-l, 7,7 ^5n-l ' " ' SfSn)- 



i=l i=l 



In the first equality, we have used the Dirac delta to integrate /i„ = 7„ ^ (7„_i/i„_]^7„_]^) • • • (71 /i^ 7^ )jn- In the 
last line, we have changed variables to gi = Jihij^ , using the translation invariancc of the Haar measure (this line 
simply is a re- writing). 

This choice of effective face weight is natural from the way spin foam models are built from BF theory. Indeed, 
the effective face weight for BF theory is a;BF(^i, • • • , ^n) = S{hi ■ ■ ■ hn)- The insertion of the group elements 7 is 
due to the simplicity constraints which break the topological nature of the theory. The way this is implemented in 
our version of the BC model is interesting because it does not change the functional form of the face weight: it is 
still formulated with a Dirac delta. A change in the face and line measures ^ would change the Dirac delta to some 
other distribution. Moreover, the form of lo allows to directly draw two conclusions. 

• Just like in the BF model, the potential divergences would come from the fact that multiplication of deltas 
may not even be defined as a distribution. More precisely, some deltas may be redundant: their arguments are 
automatically the identity of SU(2) once the other deltas are satisfiecQ. Having Dirac deltas in to makes this 
choice quite convenient to study. In particular, the fact that wbf is a delta is the reason why divergences in spin 
foams are called bubble divergences (because typically there is a redundancy for each 'bubble', i.e. independent, 
spherical, closed surface, so to speak). We will thus be able to compare the contributions of bubbles in this 
model with the BF case. 

• The expression ^ further indicates why the model is more likely to be convergent that the BF model, and 
even why faces with two vertices are the most dangerous with respect to divergences. Since SU(2) has three 
real dimensions, the Dirac delta in the effective weight ([9]) has three real components. Because the definition of 
UJ integrates them (over the conjugacy class of each hi), uj is expected to be more regular, less distributional so 
to speak, than the SU(2) delta. The more regular uj is, the more likely it is that products of uj are well-defined, 
removing divergences. Moreover, there is as many integrals in ^ as vertices around the face. Therefore, 
w certainly becomes quite regular for faces with a sufficient number of lines, and only faces with few lines 
are expected to be dangerous. We expect this feature to hold more generally in spin foam models for quantum 
gravity, as the simplicity constraints always amount to smearing the SU(2) deltas of BF theory, with one integral 
per line around each face. 



* In the third reference of [iql . it was noticed that even in the absence of redundancies, the ampUtudo may not bo finite in BF theory, 
due to singularities on the set of solutions to the constraints. This issue will be also discussed later in the paper. 



B. Analysis of the effective face weights 

1. Support and geometric interpretation 

To understand the support of a; and the content of the constraint in ([9]) , we need a bit of spherical geometrjQ. A 
spherical n-gon is a loop of n geodesic segments on the 2-sphere, with lengths in [0, n] (note that it can be degenerate 
and have self-intersections). 

Let hi, . ..,hn G SU(2) with class angles 6^. € [0,7r] defined as ^trh^ ~ cosS^. Then there exists SU(2) elements 
72 , . . . , 7n such that 

hl^2h2"f2^ ••• lnhn"/~^ =1, (13) 

if and only if there exists a spherical n-gon with lengths {9i , . . . , ^„) . The fact that the constraint implies the existence 
of a spherical polygon is proved in [32] and we will not repeat it. We will prove the reverse. The case n = 2 is trivial. 
We proceed by induction on n starting with n = 3. 

Write hk = expiOka ■ fik with Ok 7^ 0,7r, and consider a spherical triangle with lengths (01,^2, ^3)- The angles 
between the sides of the triangles are given by the spherical law of cosines, 

cos 9k — cos 9i cos 6', , , , 

cos </>.y = r-7 — r-7 ■ 14 

smt^i smt^j 

Set 72 G SU(2) such that i?(72)^2 = ^2 with 112 being any unit vector satisfying 

hi -112^ -cos 012. (15) 

The SU(2) element 72/1272^ has the same class angle as /12, but its axis is U2, i-C 72/1272^ ~ cxpi92 o ■ U2- Then the 
equation P^ for i = 1, j = 2, fc = 3 exactly reads 

tr /ii 72/1272^^ = tr /i3^, (16) 

which implies that the matrices on both sides are conjugated to each other, by, say, 73 G SU(2). In other words, there 
exist 72,73 G SU(2) such that ^1172/1272^^ — 73^3^^73^"'^- 

Completing the induction is easy. Assume there is a n-gon with spherical lengths (9\, . . . , 0„). By splitting it on its 
(n— 2)-th vertex, we get a (n— l)-gon with lengths (6*1, .. . ,6'„_2, 9q) and a spherical triangle with lengths (Sq, ^n-i, 9n). 
The induction hypothesis ensures the existence of SU(2) elements such that /ii72/i272^^ ■ • ■ 771-2^71-27^^2^0 = I where 
/iQ has class angle 9q. As for the triangle the above proof of the case n = 3 shows that there exist SU(2) elements 
such that hojn-ihn-i^n-ilnhnJn^ = ^j where ho has also class angle ^o- Therefore ho and ho are conjugated and 
we can write for the triangle /iq ^n-ihu-i^n-ilrihu^n^ = I- This completes the proof. 

Furthermore, existence of an n-gon is equivalent to the spherical polygon inequalities on {9i, . . . , 0„). This means 
that those inequalities provide the support of the effective face weight. They read 

^0,-^0,-7r(|P|-l)<O, (17) 

ieP ieP' 

for any subset PC {1, . . . , n} with \P\ odd, and P' = {1, . . . , n} \ P. 

In summary the effective face weights are supported on configurations of group elements hi, ■ ■ ■ , /i„ for which the 
corresponding class angles define a (possibly degenerate) spherical n-gon. We will find this constraint again in the 
explicit expression of the face weights. 

2. Explicit expression of the face weight 

From the bound ([5]), we see that only faces with two vertices are dangerous, suggesting that the effective face weight 
of faces with at least three vertices are quite regular. 



The authors are grateful to W. Kaminski for pointing this out. 



a. Face with two vertices, n = 2. In that case, a;(/ii,/i2) is actually well-known to be the one-dimensional delta 
constraining hi and ft.2 to lie in the same conjugacy class (i.e. to have the same rotation angle), 

uj{hi,h2)= d-y 6{hi-/h2^^^) ^y^yj{hi)xjih2) = ^—p. — r-^r V" sin fc^i sin ^6*2 . (18) 

/ ^-^ sm [71 sm C72 ^~^ 

b. Face with more than two vertices, n > 3. Since n = 2 has only a one-dimensional delta function, and there 
are more integrals over conjugacy classes for n > 3, a; does not contain any delta anymore and is regular almost 
everywhere. To make this precise, the expression of the character in terms of the class angle is inserted into ()10|) . 

^1 sin(fc6'i) •••sin(fc6'„) 
u;{hu...,h^)-)_^^^—-^——^—, 

i-— 1 

oo , (19) 

(2i)" ^^ A:"-2 ^^ sin6li---sin6l„ L ^ ^ J' 

If one group element is set to hi — ±1 (i.e. 9i = 0, tt), then w reduces to the face weight with simply one argument 
less, 71 — i- n — 1 (which is the weight for a face with one line and one vertex less). Therefore we assume that all 
h, ^ ±1. 

This allows to evaluate the sum over k for all terms independently. 



1 ^ 1 



g.fee 



^ ' J 2" ^-^ fc"-2 V 2/ 8 71-2! V27r/' ^ ' 



(2i)"^/c"-2L '^ ^ ^ J 2"^fc«-2 V 2/ 8(71-2)! 

^ ^ k=l fc=l ^ -^ 

where we have recognizeqj the Fourier expansion of the Bernoulli polynomial Bn-2, which holds for Q 6 [0, 2tt]. If Q 
is outside this interval we need to shift it back to be in [0, 27r] using that the left hand side has to be periodic in Q. 

Notice that ([20| is absolutely convergent for 7i > 4. For 7i = 3 we get X]fc°=i sin(fc0)/fc which gives the Fourier 
expansion of the sawtooth wave, i.e. a finite function which is not continous. This applies to each configuration of e^, 
appearing in ([T9)) . with &{e} = Oi + ^21=2 ^'^' mod {2tt). Thus we have 

c.(/ii,...,fe„)^ /"" ' ^ J2 e2---e^B,,J^). (21) 

8 (n — 2)1 smt^i ■ • -smt^n ^^-^ , \ 2tt / 

e2,---,e7i = ±l 

The Bernoulli polynomial Bn{x) has a monomial of highest order cc". However, the sum over the signs e^- = ± may 
lead to simplifications. 

Let us focus on the case ti = 3. There, we need the polynomial Bi{^^^^^^) = ^ + ^ — 1/2, where N € Z 
has to be chosen such that 9 + 2ttN G [0,27r]. Assume that 61,62,03 satisfy the spherical triangle inequalities, i.e. 
61+62 + 63 < 27r and 6a < 61, + 6c for any permutation of a,b,c ~ 1, 2, 3. Of the four combinations of €2, £3, three 
arguments 61 + X];=2 ^'^' ^^"^ ^^ [0, 27r]. However for £2 = £3 = —1 we obtain a negative argument. For this last 
summand we need to shift the argument back to [0, 27r] by choosing N = 1. Then in the sum over the £2, £3 all the 
linear terms in 6^ vanish and we are left with a constant 

-.(^^^ll^)-.(^^^lF^)-.(^m^)-.(^^-^..^)- ■ (-) 

Hence if the spherical triangle inequalities are satisfied, the effective face weight reduces to 

-(/.i,/.2,M==|^.^^^^J^^^.^^^. (23) 

Consistently with the support found in the Section IIIIB 11 one can check that violations of the triangle inequalities 
lead to cj = 0. Thus the face weight for a three- valent face is simply 

-(^^'^-'^-)-| sin.,sin.2sin.2 ^^^^'^"^-^ ^''^ 



We could also recognize the real and imaginary parts of the polylogarithm Lin{z) = X]fe^i ^ /k' 



where H{8i, 62,93) — 1 if the spherical triangle inequalities are satisfied and vanishing otherwise. 

For higher-valent faces, with n edges, the Heaviside like function H in 1^^ is replaced by some piecewise polynomial 
of order {n — 3) so that the support is on configurations satisfying the spherical triangle inequalities. 

As a conclusion, the face weight is a function of the class angles only. The face with two vertices is the most 
singular, and the weight becomes smoother as the number of vertices per face increases. This is expected because 
the number of group averages in the expression ([9]) for w is precisely the number of vertices. As the group averaging 
comes from the simplicity constraints, we expect this feature to also hold for the EPR/FK model. We will provide 
evidence (which will depend on the choice of certain edge weight factors) in section IIVCI 

Let us note that possible divergences of the models can be also analyzed with the help of microlocal analysis [Sj] . 
Wave front sets are a refinement (or extensions to co-tangent space) of the singular support of a function. The wave 
front sets for the effective face weights for the BC model are non-empty and indeed reflect the configurations at 



which the effective face weights are non-smooth [3J]. Thus wave front sets in spin foam models do not necessarily 
correspond to divergences but could also just signify non-smooth behaviour. In this case there is no need for regulating 
the models. 

C. Two faces glued together: the square of the face weight 

Although pretty explicit, the effective face weights are functions of the dihedral angles, which involve in their 
definition the product of two group elements: cos ^eiea = — tr h~^^he^y. Thus if we glue faces together we would have 
to disentangle these group elements again, as the /leii, and h^^v will in general be shared by different sets of faces. 

Also, even for effective face weights with valency larger than two, the sets where some hi are ±1 might lead to 
(delta-function like) singularities (this is what prevents u from being a standard function) and they may contribute, 
though having zero measure. Therefore, for practical calculations, we will rather use the expression ^ and deal with 
products of SU(2) delta functions. 

We have already emphasized that typical divergences might show up because of the products of deltas are a priori 
ill-defined, due to redundancies of the enforced constraints. Such redundancies typically arc expected when some 
faces are glued in a way that creates closed surfaces in F (boundary of 3-cells if F is the 2-skeleton of a higher 
dimensional cell complex). In the BF spin foam model, it is easy to check that the presence of 3-cells is associated 
to a redundant Dirac delta, whose argument is automatically the unit of the group whenever the deltas on the other 
faces are satisfieqj. In the following subsections, we will focus on these typical situations and argue that there is no 
divergence, expect in one case where two faces with exactly two vertices are glued together. 

The first step to study this type of situation is to look at the square of w. Geometrically, that corresponds to a 
closed surface in F made of two faces with the same boundary lines and vertices. In the BF case, cj = S and its square 
is obviously not well-defined, because the two deltas impose the same constraint twice. Here the question is therefore: 
is a; a distribution which can be squared? 

Consider that there are n lines around the two faces. There are also other faces which may share boundary lines 
with the two faces. Therefore the amplitude on F reads 

/n 
l[dK [c^(/ii,...,/i„)]'/(/ii,...,/i„), (25) 

a=l 

where / is the result of integrating all the group elements h^f in F at fixed wedge group elements hi, . . . ,h„ on the 
boundary of the two faces. / is typically a distribution, or it might itself contain divergences. However, our aim is to 
isolate the contribution of the two faces glued together and therefore we consider / as a regular function. The integral 
is on a compact manifold which implies 

/n 
l[dha [uj{hi,...,hn)]\ (26) 

a=l 

where the constant K is the maximal value of / on SU(2)". 

It is possible to calculate the integral of w^ exactly. Going through the steps of the Equation (|12p , 



6=1 



/I Ij nit I L i. 

n dha [Uj{hi, ..., /l„)] '=J l[d/3a dja J] ^gt S [{Pll^' giJlP^') {13212 '92I2P2') ' ' " l^nln'g-^i ' " ' 9T'jnPn' 

(27) 



^ The full evaluation of divergences in the BF spin foam models requires in addition to take into account the reducibility of the gauge 
symmetries and global, topological effects, see [39|. 



Re- absorbing 7^ on the right of pi gives 



a=l "^ a=l h=l 

(28) 
We then use the character expansion and exphcit integration of Wigner matrices. Doing so yields 

. n— 1 n— 1 



J]^ d/la [w(/ll, . . ■,hn)Y = Yl 3^r^ / n ^3a Xj (ff 1 " ' " ffn-l) H ^J (^a 

E ^5(^=C(2n-4) 



a=l jeN/2 J a=l o=l ,„q,^ 



,2(n-2) 
i«EN/2 "j 

This is obviously finite as soon as n > 3. One concludes that Z{r) is finite for regular enough /. 

Remark. The above result is identical to the partition function of 2d BF, i.e. 2d Yang-Mills at zero coupling, on 
a surface of genus n — 1. This is no coincidence since (P5| can actually be re- written in the typical 2d YM form. We 
first abosrb /3„ into the other Pa ■^ Pn^Pa- Then we proceed to the change of variables (ga, Pa) *~^ (kaiCta) step by 
step starting with a = 1, 

51 = 52 • • • 9n-l K^ 9n-l ■■■92^ ^ud ^1 = 52 • • • 5n-l «! g~\ ■ ■ ' 92^, 

52 = 53 • • • 9n-i k^^ g~^^ ■■■9^^ and ^2 = 53 ■ • • 5"-i "2 5^-1 • • • 53'^ 

5«-2 = 9n-i K~2 5,7-1 and /3„_2 = 5n-l "n-2 5,7-1: 

5«-i = fc,7-i and Pn-i = a,i_i. 

This recasts ((28|) in the form 

/n „ ,1-1 

]^d/la [w(/ll,...,/l,i)] =/]^(iaadfca (5([fcl , tti] [fc2 , a2] • ' ' [^n-l , an-l]) , (30) 

a=l '' a=l 

where [fc, a] = kak^^a^^ is the group commutator. The argument of the delta function in the above equation is 
recognized as a SU(2) version of the presentation of the fundamental group of the surface of genus n — 1, as expected, 
[Ml. 

IV. FINITENESS OF SINGLE BUBBLE CONTRIBUTIONS 

The calculation performed in the case of two faces glued along their boundaries generalizes to any single bubble 
of arbitrary shape. The same way we had found the partition function of 2d BF on a surface of genus the number 
of boundary lines minus 1, we will see that the calculation goes through in the case of a single bubble thanks to the 
well-known fact that lattice gauge theories are trivial in two dimensions. 

A. Single bubble contribution as two-dimensional spin foams 

Suppose we can identify in F a closed surface S of arbitrary Euler characteristic x = ^ ^ Ei + F. It is usually 
referred to in the quantum gravity literature as a (not necessarily spherical) bubble. It is such that each line is shared 
by exactly two faces. The partition function on F reads 

Z{T)=j W dhey J]^ a;(/iu;i/,- ■•>«,„ J J]^ c?/iei, J]^ ^(/i«,i/, ■ • ■ , /i!«„/)- (31) 

Integrating all the group elements h^y where e does not belong to E produces a function depending on the group 
elements on E, 

fs{{hev}e,ucs) = / ]^ d/leu J]^ a;(/itu^ ^ . . . , /l,u„j. ) . (32) 

(e.v) fglT. 
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It is a sort of Hartlc-Hawking wave-function on S but it is not clearly well-defined a priori. However, to isolate the 
contribution of S, we will consider that /s is a regular function. Because it is defined on a compact space, we assume 
it is bounded by some constant K. Therefore, Z{T) is bounded by K times the partition function where / is set to 
1. This is exactly the partition function of the Barrett-Crane model on the two-dimensional surface E. Therefore 



Z{T) < K Z(S), 



with 



•' (e,v) /CS 



«./)• 



(33) 



(34) 



{e,v) 



To evaluate the potentially divergent contribution of E C F, we will calculate ^(E) instead of Z{T). This amounts 
to ignoring the faces external to E. In the usual spin foam language, this is simply setting the spins of the external 
faces to zero. This has already been used in the literature to isolate divergent contributions. In @, it was argued 
to yield a fair evaluation of the divergences, just like in ordinary quantum field theory the loop divergences are often 
evaluated by setting the momenta of the external legs to zero. 

However, it seems to have gone un-noticed that when doing so on a single bubble, the remaining part Z{YJ) is just 
a two-dimensional version of the initial model. Thinking of spin foam models as generalized lattice gauge theories, 
and given that two-dimensional lattice gauge theories are solvable, this gives us hope to calculate Z{Y1) exactly. This 
is what we do now, in the Barrett-Crane case in the Section HVBl and for generic spin foam models in llVCI fthe BC 
case is just a particular case, but we treat them separately because the BC model is our working example in this 
paper). 



B. The Barrett-Crane model on two-dimensional surfaces 

To compute Z{T,), we use the character expansion of w, pUj) . to get 



(35) 



{i/}/CE (e.f) 



/CS 



where the sum is over all possible assignments of spins to faces, rif is the number of lines on the boundary of /, and 
/i^ is the wedge holonomy. Remember that a wedge is identified by a pair 'vertex-face', or equivalently the two lines 
along that face which meet at that vertex. We denote these two lines eu,, e^. To perform the integrals explicitly, we 
notice that they actually factorize onto vertices. 






Yldhev Ylxjfihe^yhe'^v) 



(36) 



Around each vertex with riy lines, there are also n„ faces (or rather wedges), and since each line is shared by exactly 
two wedges, we can label the lines and the wedges, say in clockwise order. At each vertex we have 



dhi--- d/i„„ Xjf, {K ^h2) Xjf^ {K ^^3) • ■ • Xj/„„ {Klhi) 



fj' 



(37) 



We have used the formula (|A6|) n^ — 1 times to integrate products of characters. As a result, all spins around v are 
must have the same value which we have denoted j. Since E is connected, the spins of all faces must be identical, so 
that the sum over all spin assignments reduces to a single sum. The summand is obtained by gathering all powers of 
dj, 

iGN/2 / V ieN/2 

On the triangulation of a surface, X) f '^/ ~ X^u '^^ ~ 2£', hence 

-2E 



(38) 



z(E) ^ y: df- 

]em/2 



(39) 
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Remark 1 — Finiteness. The result is divergent li E < x- Since x ^ 2 and we want at least E > 2, wc find that 
the only divergent case is i? = x = 2, meaning a spherical bubble with only two lines. This is the ease of two faces 
with two boundary lines glued together already seen in the section BlI CI Any other bubble is finite, 

Z(I]) = C{2E - 2x) = {-lf~-^^ ""''2 gi^^-Xr'^ • ^'^^ 

Here C, is the Riemann zeta function and Bn a Bernoulli number. 

Remark 2 — Invariance. In the BF case, the partition function on a surface of Euler characteristic x is C(~x) 
and is therefore independent of the triangulation. Here Z{T,) depends on the triangulation only through the number 
of lines and not its particular shape. It means that Z(S) is invariant under homeomorphisms of the triangulation 
which preserves the number of lines. In particular, the 2-2 Pachner move does so and therefore leaves Z(T,) invariant. 

Remark 3 — 4-2 Pachner move. The section IIII CI is obviously a particular case of this section, with x — "^ 
and E = n the number of boundary lines of the two faces. A more interesting case for quantum gravity is the 
4-2 Pachner move. As a Pachner move, it is a change of triangulation which preserves the topology of T. In the 
dual of F, one changes a configuration of two 4-simplices which share a common tetrahedron and thus have eight 
boundary tetrahedra, with four 4-simplices each contributing to two boundary tetrahcdra. The four 4-simplices are 
glued together in a specific pattern such that two 4-simplices share exactly one tetrahedron. In the 2-complex F, 
we find 4 vertices completely connected by 6 lines which form 4 triangular faces. This pattern corresponds to the 
boundary of a tetrahedron and this is precisely the surface S, with % = 2, i? = 6. Thus the contribution to the 
partition function from the bubble in the 4-2 Pachner move configuration is finite. 

C. Single bubble contributions for general models 

Spin foam models in 2D reduce to 2D (standard) lattice gauge theories [7|] , implying that the above calculation for 
a single bubble can be extended to generic models. A wide class of models (including BC with generic choices of edge 
and face weights and the EPRL/FK model) Q is given by 

Z= f l[dh,f'[[C{{Kf}f^,)'[[w{h,,f-K^j) . (41) 

^ [ej) e / 

Here we integrate over group G elements associated to edge- face pairs (e/). For every edge we have an edge weight 
C which depends on the group elements /le/, for which e is an edge in the boundary of /. For each face we have a 
('bare') face weight w , which is a class function and evaluated on the holonomy around the face. 

For a 2D surface, there are always two faces adjacent to a given edge. The edge weights C have also to satisfy a 
certain invariance property, which means that for the 2D case we can expand C as follows into irreducible unitary 
representations p of G 

C(/ii,/i2)=^C'pdim(p)xp(/ii/i2"') . (42) 

p 

The face weights are expanded as 

z«(M=^(ipdim(p)xp(M . (43) 

p 

Using these expansions in (|4ip for the 2D case one notices that the sums over the representation labels p reduce 
to one sum as the group integrations impose Kronecker deltas between representation labels. Taking care of all the 
dimension factors that come from the expansions (|42l43p . the group inner product between representation matrix 
elements, as well as from the contractions of Kronecker deltas (which gives traces around vertices) we obtain 

Z2D = ^(Cp)^((ip)^(dimp)^-^+^ = ^((7p)^(cip)^(dimp)'^ . (44) 

p p 

For (standard) lattice gauge theory we have Cp = 1, and we recover the corresponding 2D partition function. For 
the BC model we have G = SU{2) x SU{2) and with our choice of edge and face weights 

'^P = l : C'p = Jp.(,-,)(j-j (45) 
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so that we recover ((39)) . However we see that changing edge weights, for instance introducing a factor d3 per edge, so 

that Cp = Spijj\ would lead to a triangulation invariant but divergent result for a spherical bubble. 

For the EPRL model (G = SU{2) x SU{2)) with Barbcro-Imniirzi parameter 7, leaving the face weights Qp free 
for the moment, we have for the edges [7| 

~ Y^ dim(j) /1+^ ■ |l"^l ■^^ /J ^ ^^2 /.rn 



d\in{—j^j) dim( 



2 



J) 



Here (ie(/9) is an edge weight factor which is left undetermined. We obtain a non-vanishing coefficient only if the 
SU{2) X SU(2) representation p = (/, j") is of the form {^j, ^-^j). 

We have two free functions, the face weights uip and the edge weight factors de{p). The choice of these factors 
will heavily influence the convergence properties of the model. Different requirements have been proposed to fix 
these weights [12h14| |. With the simple arguments put forward here we can comment on how these requirements will 
influence the convergence of single spherical bubbles. 

(i) We can require that the model is invariant under edge subdivisions (here for edges which are shared by only two 
faces) and face subdivisions. To achieve this we choose Qp = 1 and rfe(p) such that Cp equal to one or is vanishing 
(if p is not admissible). This will give a triangulation invariant 2D model. In this case spherical bubbles will diverge 
(assuming that there are infinitely many admissible representations). 

(ii) A weaker requirement is invariance under subdivision of faccs|f| Dividing one face into two we raise the number 
of faces and the number of edges by one. Thus invariance requires uip = (Cp)^^ for admissible representations p. The 
convergence then depends on the difference between the number of faces and edges. If we consider the square of an 
effective face weight, it forms a spherical bubble with two faces. As long as Cp scales with some negative power of 
dimp we obtain a more convergent result with growing number of edges, where the specifics again depend on the edge 
weight factor. 

If LOp = (Cp)~^ we have a divergent partition functions for all spheres where the number of edges and faces are 
equal to each other. This includes the bubbles that appear in the 4-dipole configurations. 

(iii) One can also adjust the face weight and edge weight factors to obtain convergent results for specific families of 
bubbles. 

Thus wc sec that we can get an estimate (as wc ignore the contribution of faces which connect to but arc not part 
of the bubble) on the behaviour of bubbles and the behaviour of effective face weights by quite simple methods. This 
allows to choose the edge weight factors and face weights according to the divergent or convergent behaviour one 
wants to achieve. We have seen however that requiring triangulation independence in 2D (i.e. invariance under face 
and edge subdivisions) comes at the cost of divergent spherical bubbles and distributional effective face weights for 
arbitrary number of edges. 

V. MULTIPLE BUBBLES: THE 4-DIPOLE 

Next we will discuss a configuration with multiple bubbles, known as the 4-dipole. It in particular it arises in 
discussions of group field theories as configurations which include such dipoles are the most divergent ones [40| . For 
the BC model the 4-dipole could be divergent due to the appearance of two-valent faces, which glue to (multiple) 
spherical bubbles. These bubbles have three faces and three edges - thus the single bubble contribution as discussed 
in section IIV Al converge. The assumption of this section was to ignore faces external to the bubble. Thus this is 
also a test whether this assumption holds in this case. Indeed, we will find that divergences do not occur - at least 
no divergences due to redundancies of delta functions. (There are however singularities on a measure zero set and we 
leave the integrability of these singularities open.) 

An expansion in spin variables for this case is not sufficient to determine convergence. We will therefore switch to 
the group representation and basically show that no redundancies arise if we solve for all the delta-functions appearing 
in the partition function for this configuration (the method is detailed in p^). 



Here we mean that a face is subdivided by a two-valent edge which goes between two already existing vertices, i.e. we do not create new 
vertices. We should point out that there exist other notions of face subdivisions which create new vertices and hence a larger number 
of additional edges. 
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Figure 1. The vertices A, B are dual to 4-siinplices, which are glued together by four tetrahedra here represented as the lines 
1,2,3,4. The lines 0^4,05 stand for the two boundary tetrahedra. 

Let us describe the 4-dipole configuration. Wc consider a piece of triangulation with two 4-simplices which are 
glued together along four of their five boundary tetrahedra. These four internal tetrahedra have six triangles and four 
edges in total, which are thus shared by the two 4-siniplices. The boundary of the gluing consists of two tetrahedra. 
Notice that their triangles are shared with internal tetrahedra (in a 4-simplex, a triangle is shared by exactly two 
tetrahedra). Therefore, each triangle of a boundary tetrahedron is identified with a triangle of the other boundary 
tetrahedron. This means that the two 4-sinipliccs actually share all their triangles. 

To write the spin foam amplitude, it is easier to work in the dual of the triangulation, depicted in the Figure [TJ It 
is called the 4-dipole because it has two vertices (of degree 5) connected by four lines, with one external line hanging 
out of each vertex. The vertices are labeled A, B and the internal lines 1, 2, 3, 4. The two external lines are denoted 
OA and OB. The lines 1,2,3,4 create six internal faces, labeled by the pairs of lines, (ij), 1 < « < j < 4, since each 
face goes along two lines only. The four edges shared by the two 4-simplices correspond to four bubbles in F, whose 
boundary are the faces (ij), (jk), (ki), for 1 < i < j < fc < 4, glued two by two. The surfaces of the bubbles are 
therefore spherical. 

The triangles of the boundary tetrahedra are dual to external faces. As the two 4-simpliccs share these triangles, 
the external faces actually go along both A and B. There are four external faces, all going along the line OA, then 
choosing an internal line i = 1,2,3,4, and then all going along the line OB. Note that these faces are broken (they 
are not closed, because the dual triangles are on the boundary of the gluing). 

In the spin representation, the amplitude for the 4-dipole has two lOj-symbol, one associated to the vertex A and 
the other to B, and they depend on the spins associated to the faces. Since the corresponding 4-simplices share all 
their triangles, it means that the spins on the faces are all common to the two lOj-symbols (there are actually only 
ten faces in the 4-dipole). We thus get the square of a lOj-symbol. Moreover, the spins of the six internal faces must 
be summed, while the spins on the external (broken) faces are fixed to ji, i = 1, 2, 3, 4. The partition function for this 
piece of triangulation is a function of the four external spins. 



^4-dipolo(jl,J2,J3,j4) =: ^ 



l<i<j<4 




(47) 



There is no dj factors because each face only goes along two lines. 

To estimate the potential divergence of these sums, it has been proposed in [9| to set the spins of the four external 
faces to zero. With ji ~ 0, each lOj reduces to the square of a Wigner 6j-symbol (this can be seen in ((!]) which directly 
reduces to the group integral formulation of the square of the 6j-synibol) . Therefore the amplitude is 



^4-dipolc(0, 0,0,0) 



E 

J12 ,jl3 J'l4 
J23 j'24 J34 



Jl2 jl3 .714 1 
J34 J24 .723/ 



(48) 



We consider the large spin behavior of the summand, when all spins are homogeneously scaled by A 3> 1. The 
Ponzano-Regge asymptotics of the 6j-symbol states that {6j} ^ 1/A^/^. The summand thus behaves like 1/A^ and 
as there are exactly six sums to perform, it is not possible to conclude about the convergence/divergence of the 
amplitude. This shows that the choice of measure in our version of the BC model requires a more subtle analysis. 
However, if gauge symmetries, or redundancies of delta functions in the group integral formulation, were present, a 
positive exponent of A would certainly be expected. It means that the simple power-counting argument still suggests 
the absence of redundancies and gauge symmetries. This is what we will show below. 
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This is done by moving first to the group integral formulation with effective face weights. The faces (ij), for 
1 < i < j ^ 4, have effective face weights 

rf7u S [h^A hj} jij hjB h-^ -i-^y (49) 

As the external faces are not closed, we cannot use effective face weights. But using the definition ^ for the two 
lOj-symbols, it is easy to see that 

» 4 4 

Zi-dipo\c{ii,J2,h,k)= j Wd-HjY\_dhiAdh,B W^yh^Ah^l lijhjBh~B I'^^j Wxni^oAKl) Xni^^BKB) ■ (50) 

•^ i<j 0=1 i<j i=l 

In other words, in the character expansion ([TU| of the external faces, only the characters on the wedges of A and 
B appear. Since there are no delta functions on the external faces (the reason being that we work at fixed external 
spins), the bubble divergences can only come from the product of delta functions on the internal faces. Therefore, we 
simplify the analysis by removing the contribution of the external faces, as in the Section lIV A[ which here amounts 
to simply ignore the oscillations of the characters Xji by setting ji = 0, like in Q. 
The partition function we want to evaluate is 

/■ ^ 

•Z'4-dipoio(0, 0, 0, 0) = / J]^ d'yij J]^ dh^Adh^B J]^ ^f /itA hj^ jij hjB Kb Itj^j ■ (51) 

By a redefinition of the elements hiB <— hiBh^g and hiA ^ hiAh^A' * ~ 2,3,4, the elements hiB,hiA are trivially 
eliminated. Then the faces (12), (13), (14) arc used to integrate hiA, i = 2,3,4, which are constrained to be 

hiA = In hiB lii^ , 1 = 2,3,4. (52) 

Only three delta functions remain, which impose constraints between the elements jij and h2B,h3B,h4B- Dropping 
the subscript B, 

^4-dipoio(0,0, 0,0) = / d/i2rf/i3rf^4 JJ^Tu i5(7i2/i27iV 713^3^^1^723/13^2^^23^) 

Kli2h2li2 liiK'^lil 724^-4/12^^2^4^ )'^(7i3/i37r3^ luK^lu l^h^h^'^-i^l). (53) 
712 can be completely absorbed into a re-definition of the other 7^ <— 7^^ 7^ . We obtain 

^4-dipoio(0,0, 0,0) = / dh2dh3dh4d'yi3d'yi4d-f23d-f24d'y34 <5(/i2 7i3/i:rSiV 723/13/1^^23^) 

S{h2 714/ir Vil^ 724/14/1^^24^) Hlishsirs luK^lii 734/14/1:^^34^ )• (54) 

Notice that when a group element appears in one delta, it also appears with its inverse in the same delta. It means 
that it is now necessary to study the remaining constraints in depth. Let us look at the constraint on the face (23), 
appearing in the first line of the above Equation. It reads 

/12 llzK^lli = 723/l2/l3"V2V- (55) 

We write hi = cos 9i + isinOihi ■ a and take the trace (in the fundamental representation) on both sides. Using the 
fact that tv{hg) = cos0hcos9g — sin 6 h sin 9 g{'hh ■ ng), it becomes when 02,^3 7^ 0,7r, 

^2 • -R(7i3)«3 = "2 • ^13, (56) 

where i?(7i3) is the matrix of 713 in the 3-dimensional, vector representation. The solution of this equations are easily 
found for arbitrary n2,n3, 

__ -i4>i3h2-ff iearccos(n2-n3)("2Xn3)-(? ^iOi^fis-a (fi7) 

The angles (pia, 613 are totally free and parametrize arbitrary rotations around n2, i^3- Notice that e = 0, 1 is a discrete 
ambiguity due to the fact that 71,3 can be rotated in the plane spanned by n2, n^ by an angle which is twice the angle 
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between n2 and fi^ without changing the scalar product n2 ■ n^. Since it is a discrete ambiguity, we can restrict 
attention to the case e = 0. Inserting this solution in the initial constraint, we find that h2h'^ has to commute 
with e"^""^'^723. Two commuting SU(2) elements must lie in the same U(l) sub-group generated by their common 
rotation axis. As a result, 

723 = e-*'^""2-^ ^ie2zA23-s^ (58) 

where 7123 is the rotation axis of h2h'^ ■ The constraint ([55]) thus admits solutions for arbitrary generic h2,h^, where 
among the six real degrees of freedom of 713,723, only 013,^13,^23 are left undetermined. It means that three real 
parameters have been fixed by the three real constraints. 

The same reasoning applies to the constraint /i2 714/^4 7]^ = 724^2^4 724 • It leads to solutions for arbitrary ft,2, ^4 
and 

up to some discrete ambiguity, where the angles (/>i4, 0i4, 024 are free. The final step is to insert these solutions into 
the last set of three real constraints, 713/137^^ 714^4 7i4 734/14^3^ 734 = I, and see whether none of them are trivially 
satisfied. The following projection of the constraint 

tr 7i3/i37r3^ 714^4 Sr4^ = tr h^ h^^ , (60) 

gives R{e-'"^^^^'-^)n3 ■ i?(e-*'^""^'')fi4 = m ■ m, which is solved by 

gi(0i3-0l4)n2-i? _ g-#34n3-'? ^irp3ifi-4-3 /Q^\ 

for some angles (/)34, tpsi, up to some discrete ambiguity. For generic /12, /13, /14, the rotation axes n2,fi3, fii are linearly 
independent. Therefore, there is no non-trivial solution in the neighbourhood of the trivial solution 0i3 = (/)i4, ^34 = 
^"34 = 0. It means that (150]) indeed removes one degree of freedom by setting ^13 ~ (f)i4. Finally, the constraint 
imposes that h^h^^ commutes with e"'^^^"^'°'734. This fixes the two real degrees of freedom of the rotation axis of 
gi0i3ra2o-^g^ and leaves one angle, denoted ^34, free, 

734 = e-*'^""^-^e**3^"='*-^, (62) 

where 7134 is the rotation axis of /13/14 . The last set of constraints has thus eliminated three real parameters, meaning 
that we did not meet any redundancies while solving the constraints. 

The parameters left undetermined are ft,2, ^13, h^, ^13, ^13, ^23, ^14, ^24, ^34, and the set of solutions is 

^ = {h2, hs, hi, 713 = fc2 e*^-"-^ 714 = fc2 e*^""-^ 723 = fc2 e*«-"--^ 724 = fc2 e*''-"--^ 734 = fc2 e^'^-"-'^; 

with /i2, hs, /14 arbitrary and ^2 = e"'"''"^-'^} . (63) 

This is a fifteen-dimensional space, while there were 8 SU(2) elements to integrate in (|54p . i.e. 8 x 3 = 24 real variables. 

Therefore, the constraints have put restrictions on 24 — 15 = 9 variables, corresponding to the number of constraints. 

This means that there are no divergences coming from the product of delta functions. Defining the map H : 

SU(2)8 ^ SU(2)3 

H{h2, hs, /14, 713, 714, 723, 724, 734) 

= (/l2 713/I3 ^7lY 723/13/1^^72^' /12 luK^J^i^ J24h4h^^-f24, 713/1371V 714/11^1!^ 734/14/13' ^73^4^) : (64) 

the above analysis reveals that for a generic solution (p Q T, dimkeriii?^ = 15 and ikdH^ — 9. Therefore the tangent 
space at cj) decomposes as T^ SU(2)* = T^J^ © N^J^, where the normal space N^T is the ortho-complement of the 
tangent space to the space of solutions. The restriction dH^^^^jr to the normal space is an invertible map from N^T 
to T^ijj) SU(2)'^. The integral over the normal directions correspond to the parameters which are fixed by constraints. 
The partition function becomes 

^4-dipoic(0, 0, 0, 0) = / dcj, ^ -. (65) 

Jjr \<iet dH^lN^jr\ 

For generic solution (j), this determinant is non-vanishing, but it may happen that it actually vanishes on a subset 
of measure zero in J-'. Such singularities are well-known to arise in BF theory [18| . In two-dimensional BF, there 



are such singularities on the moduli space of flat connections |38j , but it can be shown that these are integrable (as 
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expected since the partition function can be exactly calculated using the character expansion as C{2g — 2) on a surface 
of genus g). However, there is absolutely no generic result on such singularities beyond the two-dimensional case, 
meaning that one has to deal with them case by case. In the present case, this is quite complicated because of the 
number of variables involved. We will not discuss further the possibility that the integral is divergent due to such 
singularities. Instead, we conclude this Section by emphasizing the fact that there are no redundancies in the initial 
product of delta functions, removing the expected source of divergences in spin foams. 

VI. GENERIC EVALUATION OF THE BC PARTITION FUNCTION 

The technique used in the last section was developed in [l8| to integrate over the representation variety of finitely 
presented fundamental groups. It can also be applied to the BC model on generic 2-complexes. 
Combining (J8]) with ([9|), we see that the partition function writes as an integral over 

A= \^A= {{hev}eeruvero, {lvf}ve^ojer2)] , (66) 

subjected to the constraint 

'- J /ei 2 

where ui, . . . , w„ and ei, . . . , e„ are the vertices and lines around the boundary of each face / G r2. -ff is a map from 
^to SU(2)|i^2l and 

Z{r)= [ dA 5{H{A)). (68) 

J A 

6{H) is the 3|r2|-dimensional delta function over the target space SU(2)l'"2l, 5{H) — Y\fS{Hf). Thinking of the group 
elements as holonomies of a gauge field, and of Hf as the corresponding Wilson loops, the integral corresponds to a 
lattice gauge theory on a two-complex at zero coupling. The integral localizes on the set J' = H^^{T) (which is the 
set of flat connections in the lattice gauge theory interpretation) . 

Therefore, we have to solve the constraints, meaning that we need to find among the real degrees of freedom of 
A G A those which are free and parametrize the set of solutions J-, and those which arc functions of the free parameters 
as determined by the constraints. Because {Hf = I)^gr2 is a set of polynomial equations on SU(2), the set of solutions 
J^ is a real algebraic variety whose dimension is the number of free parameters. For generic solutions G J-", this 
turns out to coincide with the dimension of the kernel of dH^. This is because 

T^J^ = keTdH^, (69) 

as expected. 

Now we can understand the potential divergences coming from the product of delta functions. Notice that the 
rank of dHcj, is ikdH^ ~ 3|r2| — kerdiJ^ and corresponds to the number of directions spanned by dH^ in the target 
space SU(2)l'"2l. If rkdi?^ is strictly less than dimSU(2)''"=l = 3|r2| for generic <^, then it means that some directions 
are not explored whatever the variations around are. Therefore the components of the delta functions along these 
directions are trivially satisfied and the amplitude is divergent. 

It is actually possible to describe the divergence rate, following [l^ [33 . If the delta functions are regularized with 
a thin width 1/A (using a heat kernel for instance), then the divergence degree is A'^l'"^'"''^'^^*. In the case of BF 
theory on a 2-complex F, it is possible to relate this divergence rate to the topology of F (and even to the spacetime 
topology if F is the 2-skeleton of a cell decomposition of a four-dimensional manifold). However, in the Barrett-Crane, 
we have not found a simple topological interpretation of the divergence rate 3|F2| — rkdiJ^. 

When rk dH^ is exactly the dimension of the target space for generic (/>, one can conclude that there is no divergence 
coming from the product of the delta functions, i.e. all the constraints are independent. This was the case for the 
4-dipole. Moreover, the integral 2'(F) rewrites as an integral over J^(F), and for each £ -^(r), an integral over the 
directions orthogonal to T^T^V)^ denoted N^J- for 'normal space'. The directions of the normal space are those along 
which the constraints fix the variations around to vanish. Therefore, the integral becomes 

Z(F) = / d^l da 5{dH^{a)) = f dc^ — — i , . (70) 

dHii,\pf jr is the restriction of dH^ to the normal space (intuitively, its kernel has been removed). 
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VII. SOME GAUGE SYMMETRIES AND RECURSION RELATIONS ON THE lOJ-SYMBOL 

A. Existence of gauge symmetries at certain solutions 

Gauge symmetries around a solution G J^ correspond to directions which are not spanned by dH^. To see their 
action, let us re-write the 5{dH^{a)) of ([7D|) as 



5{dH^{a)):^ db cxpi{b,dH^{a)), (71) 

J5u(2)ir2 

where (•, •) = J2f{'j')f i^ the sum of the invariant inner product over tlic different copies of su(2). The variable 
b = {bf}f er2 is a- Lagrange multiplier imposing the constraint. A gauge symmetry is a i/i-dependcnt, non-zero 
variation 6^b which leaves the action {b,dH^{a)) invariant, 

{S^b,dH4a))^0, (72) 

for any variation a £ T^A. Typically, we want a gauge symmetry to involve non-trivial 5^bf 7^ on several faceo 

When the constraints are independent, diJ^ has maximal rank and there is generically no gauge symmetry. We 
have seen this is the case for the 4-dipole, and single bubble contributions are finite so they have no gauge symmetry 
either. This is expected in theory of discretized gravity such as spin foams, because going on the lattice breaks 
diffcomorphism invariance [26| . There is however a special class of triangulations which only admit flat solutions and 
for which one would expect diffeomorphism symmetry also in the discrete case [4l|. Physically, the gauge symmetry 
in the flat space case means that vertices of the triangulation can be moved around without changing the physics. 
Such a vertex translation symmetry arises for instance for the 5-1 move configuration (arising from a subdivision of a 
4~-simplex into 5 simplices), where the inner vertex can be moved in four directions without changing the fiatness of 
the configuration. While these symmetries have been observed and (canonically) analyzed in details in Regge calculus 
[231 a-iid shown to lead to Dirac's hypersurface deformation algebras [23, the (quantum) spin foam case remains 
mostly unexplored (with the exception of 3D gravity and BF theory in higher dimensions [13, 123, UM ■ Assuming such 
symmetries would exist, we could derive Hamiltonian constraint operators for the boundary wave functions defined 
by the spin foam transition amplitudes. This would lead to a canonical theory describing the amplitudes defined 
by spin foams and thus to a connection between canonical LQG and spin foams. The case of BF theory, where the 
symmetries are not broken, has been discussed in |22| - [24| . For the BC model we found only one configuration where 
such a symmetry occurred - the case of two two-valent faces glued to each other. Indeed this constitutes the only 
bubble divergence we found. Hence we do not expect a divergence (and hence no full symmetry) for the 5-1 move 
configuration. 

Additionally to the special class of triangulations which only support flat solutions, such vertex translation sym- 
metries might arise around special (i.e. flat) solutions in more general triangulations. The Hessian evaluated on such 



solutions will have zero modes corresponding to infinitesimal vertex translation symmetry [421 ] . Again such symmetries 
have not been discussed in the spin foam case yet. Here we ask whether a possible similar phenomenon exist for spin 
foams, i.e. symmetries which occur only at special solutions. 

Indeed we will show that some solutions (j) € T to the constraint have gauge symmetries. We have already mentioned 
that even when the constraints are independent for generic solutions, there may be some solutions (a set of measure 
zero in J-) which are singular because detdH^ij^jr = 0, meaning that "kcvdH^ becomes larger than T^J- and not all 
directions of the target space SU(2)I'"2' are spanned. Some of these singularities may be interpreted as the appearance 
of gauge symmetries (j72p (but not all singularities correspond to gauge symmetries, since the phenomenon of footnote 
|9]can also occur in singularities). 

The condition (|72p has to hold for any variation, in particular when applied to a variation ^^^y of a group element 

(5^6/, Ad(7^,i//leit-i/i;:2l,i7j;/ • ■ ■lv,-lfhe,-lviK^^lv^_^f) 1 - ^^ilv^f^e^fK^^flv^) ^v,f) = 0" (^4) 

There are two obvious situations where this is true: 



In BF theory on the torus, 



: / dxdy 5{xyx-^y-^), (73) 



/SU(2)2 

the constraint forces x and y to He in the same U(l) subgroup, say generated by tr ■ n. Then it is quite easy to see that the hnearized 
constraint d{xyx~^y~^) never spans the direction ct ■ n G 0u(2). However, this is clearly not a phenomenon wc want to call gauge 
symmetry. 
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• when all the group elements h^v are the same around each vertex hg^ — /i„, because the operator 1 — 
Ad{"fTj.fhf..fh~ flyf) is then identically zero, for any elements jyf, (below we will need in addition that 
the elements 7^/ only depend on the face), 

• when all group elements 7^/, hev lie in the same U(l) sub-group of SU(2), say generated by c? • n, and S^b are 
variations in this direction. This is because [1 — Ad(e*""'^)](n • (?) = 0. 

In both cases, we can exhibit gauge symmetries, when F is a d-dimensional cell complex, with d > 3, using cellular 
homology. The chain spaces Ci(r) are real vector spaces, Ci ~ IRI'^'I (hence we identify chains and co-chains), with 
boundary operators dt and co-boundary operators (5*, 

5° 51 s^ s"-^ 

0^ Co(r) ^ Ci(r) ^ C2{T) ^ ••• ^ Cd(r) ^o, (75) 

di 82 ds dd 

The boundary and co-boundary operators satisfy di^i o di = and (5' o J'^^ = and they are dual to each other. 
The operator S^ sends lines to faces: ii v = {we}eeri then S^{v) = {<5^(w)|/}/gr2 with 

'5'(«)|/ = E^-/«-' (76) 

eC/ 

where eg/ = ± denotes the relative orientation between e and /. 
In the case the elements h^v only depend on the vertices, 

hev=K, and 7„j = 7^, (77) 

on the faces, the differential of H reduces to 



Ad(7/)Eee/ (c«/le.(e)/l4e) -«'^et(e)/iet(e))■ 



(78) 



Here s(e),i(e) denote the source and target vertices of the line e. Therefore, S^ and dH are simply related. Denote 
Qh =■ dhh^^ : Tu SU(2) -^ su(2) the Maurer-Cartan form which maps the tangent space at h to the Lie algebra. It 
becomes 

rfi?{/i„=;i„,7„^=7^}|/({ae«}) = -(5}(8) Ad(7/)eft„(aes(e) -a^t(e)), (79) 

for any tangent vectors Oeu G T^^ SU(2). The adjoint action by 7/ on each face can be absorbed into a re-definition 
of the Lagrange multipliers bf <— Ad(77^)6/, so that dH is basically the cellular co-boundary operator 5^ . We notice 
that 

{b + di®iA,^(2){c),dH^{a)) = {b,dH^{a)) + (c, (<5^ ® id,u(2)) o diJ^(a)), (80) 

using (93 ® idsu(2)(c), x) = (c, S^ ® id5u(2)(a::)). Thanks to the identity S^ o S^ = 0, we see that 

6^5 + a3«)id^^2)(c), (81) 

for any c G C'3(r) C>$su(2) is a gauge transformation. This is the same gauge symmetry as in the topological BF theory 
with structure group R'^. The reason is that the solution we are looking at is up to local rotations (at the vertices) 
equivalent to the trivial solution where all group elements are the identity, and in the neighborhood of the identity, 
SU(2) looks like R^. 

In the case where all group elements are generated by a single direction, a ■ n, the gauge symmetry is the same as 
in a U(l) BF theory. Let us parametrize the group elements as 

The constraint Hf{A) — I then reduces to a U(l) constraint, 

E Oesie) - Oetie) = mod (2^), (83) 

eCf 
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making the contact with U(l) BF theory obvious. The privileged direction a ■ h induces a natural splitting of 
su(2) = u(l)Ti ® Q± where u(l)^ = spanja • h} and q± is its ortho-complement. This also gives a natural basis in 

rft„SU(2), 

T^,„ SU(2) = Rae„ © 0e„± , (84) 

where Qev± is spanned by the derivatives with respect to the two components of the rotation axis n^v of h^v evaluated 
at fieu = h. A similar decomposition T^^^ SU(2) = ]R9q^,^ © 9vf± holds. 

Some straightforward algebra shows that dH sends dg^^ to u{l)n, and Qev± as well as 2vf± to g±, and that it 
vanishes on da^f ■ Moreover, the restriction of dH to the sub-spaces dg^^ basically reduces to the cellular co-boundary 
operator S^, 

dHf{{Xevde^^}) = -SJ{{Xes{e) " Xet{e)}) ® S ■ fl. (85) 

Therefore, the action is left invariant by the transformation 

&^6 + 93®id.u(2)(c), VceC3(r)®u(l)ft. (86) 

Since these gauge symmetries rely on cellular homology, they are reducible as soon as d > 4. Indeed, the gauge 
parameters c are not independent. If two of them differ by di{y) for y e C'4(r) ® su(2) in the first case and 
y € C'4(r) © u(l)ft in the second case, then they induce exactly the same gauge transformation (because 930^4 = 0). 
This reducibility is well-known in BF theory [l7[. 

B. Recursion relations on the lOj-symbol 

1. Using the 4-dipole 

We consider the 4-dipole configuration as in the Section |V] where we wrote the partition function ([50]) with fixed 
spins on the external faces. However, in this partition function, not all solutions are of the form ((77)) . Further, we 
have seen that there is no gauge symmetry for generic solutions since there is no redundancies in the constraints. 
Therefore, instead of the partition function Z4_dipoio in (|50p . we will consider the following quantity, 

/4-dipolo(jlB,J2B, J3B,J4B,J23, J24,J34) = / J^d^J/ld/ijB J| ^7^ ^ (^(^OA/liyi^) Xj.b (^lOB^iS ) 

•^ i=0 l<i<J<4 i=l 

Yl S{hiAhJ^-fijhjBh^Bli/) n XhA'^i^'^jh^jhjBK^-f:r.^), (87) 

j=2,3,4 2<i<j<4 

and proceed to evaluate /4-dipoic in two different ways to get recursion relations on the lOj-symbol, which can be 
interpreted as a (constraint) equation on the vertex amplitude. 

But first let us point out the difference between Z^4_dipoie and /4-dipoic- First, note that the integration variables 
are the same, only the integrands and the boundary variables differ. We have changed by hand the characters of the 
external faces at vertex A in (|5(H like 



444 
Yl^nihoAh-2) ^ l[^{hoAK2)= Y. li'^nA Xj.AhoAh-2). (88) 

i=l i=l jlA,J2A i=l 

33A,J4A 

This allows to satisfy the special condition ([77)) on the elements hev at the vertex A. Notice that the delta functions 
have an expansion onto characters similar to the initial characters of ()50p . However, in Z^^dipoic, the spins of the 
characters Xjii^oAh^A) ^^ the vertex A and of the characters XjtihoBh^B) at the vertex B are the same, because 
they correspond to the same (external) faces going along both A and B. When putting in /4_dipoiG some additional 
constraints, the equality jiA = jiB = ji is broken because more modes are necessary to enforce the condition ()77[) at 
vertex A. 

The dipole possesses six internal faces. Combining the constraint (|55)) with the effective face weights of the faces 
(li), for i = 2,3, 4, imposes in turn that hiB = hjB, for any 1 < i < J < 4. Indeed, 

S{hoAh^l) S{hoAh~2) / d-fu S{hiAh~2'lijhiBh^B^^^) = '^(^oA/ii'i) H^-oAh'^) SihiBh^s)- (^9) 
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As for the three remaining faces (23), (24), (34), the delta functions of their effective face weights is automaticaUy sat- 



isfied thanks to (IMl) . i.e. /i,;4/i„ l^nhjBh^g^.^,^ = I (for any 7^). Therefore, these delta functions become redundant. 



iA'I'jA'yij 



confirming in this case the existence of gauge symmetries. To avoid the divergences associated to these redundancies, 
observe that we have not included the effective face weights of the faces (23), (24), (34) in /4_dipoio- Instead, we have 
only picked up one mode of their character expansion (the last line of products in ((87)) ). Due to (|89p . these characters 
simply evaluates to the dimension of their representation, 



Y[S{hoAhJ) 6{h,B Ihe) H ^J^ji'^^Ahjl-frjhjBh^il,/) ==YlS{hoAhJ) S{h,Bhj^B) Yl 



2<i<j<4 



i=l 



2<i<j<4 



(90) 



This product of dimensions is the sole dependence of /4-dipoic in the spins j23,.724, ^34. If these spins were summed 
(with measure dj-^ ) to form the effective face weight as in .^4_dipoioi we would get X^i "^^ which is obviously divergent 
(it is the formal expansion of (5(1)). The fact that the dependence of /4-dipoio on jij (2 < z < j < 4) is just dj.- is the 
signature of the gauge symmetry, similarly to the case of spherical bubbles in BF theory |23 | 

Therefore, the only non-trivial contribution to /4-dipoic is the product of characters on the external faces at the 
vertex B. With the change of variable h — hoBh^B (this quantity is independent of i = 1, 2, 3, 4), we finally get 



'4— dipolc 



^J23 %24 "J34 



dh Ylxj.nih), 



(91) 



The second way to evaluate /4-dipoio is through a character expansion of all the delta functions, and integrating 
the variables 7^^ . That leads to 



^4— dipolc 



E 



JlA,]2A,33A,JiA 
il2 J13 J'l4 



^J23 ^J24 "•. 



i=0 l<i<j<4 



/- 4 _ 4 

l[[dh,A n X],A^^Ahj})Y[x].AihoAhJ 

■ 4 ., 

jWdhiB Jl Xj,,{htBh^B)W_X],B{hoBh^B) 



i=0 l<i<j<A 



(92) 



The two quantities into square brackets are lOj-symbols, according to the definition ([TJ. Equating this formula with 
(|9T]) leads to 



E 



n-2-jii-ni 



E 



n^^- 



jlA:i2A V.i=l 
J3A ,UA 





[dj23 d32i dj^^f dh YIx^bW. (93) 



This is a new sum rule for lOj-symbols. The summand is not the square of the lOj-symbol as in the spin foam 
representation ([M)) of ^4_dipoie because of the constraints (|55)) which implies jiA y^ jiB- Moreover, the spins jiA are 
summed here while the boundary spins ji in Z4_dipoie are fixed. The final difference is that the spins J23t J24-, jsi are 
arbitrary but fixed, instead of being summed, to avoid the divergences due to the gauge symmetries. 

Our formula can be specialized to specific values of the spins free spins jiB^i = li 2,3,4, j23,j24:j34- For instance, 
setting jiB = 0, the lOj-symbol in the second bracket collapses to a squared 6j-symbol, and the integral on the right 
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hand side is trivialized as xo(^) ~ 1, 



E(jl2 .713 juV 
. . . 1J34 J24 J23j 
J12J13J14 



E 



n^^- 



J1A,J2A Li = l 
J3A,jiA 




[dj23 Ctj24 '^J34j 



(94) 



Another interesting way to use our main formula (|93p is to sum over one of the spins jiB, say Jib, with measure 



dj-^g . Then the integral on the right hand side simplifies, 



Y^ dj,g j dhW XnB {h) ^ dh S{h) Xhn W X^b W Xmb W == d^^^ dj.,^ dj^^ . 



(95) 



Therefore, 









>; 

Jl2,jl3,jl4, 


>; 

jlA,J2A 
J3A ,J4A 


[[d,,A 




E^. 




— [«J23 "'J2i d 



'J34J 



"j2B dj3B dj4,B ■ (90 j 



One can then further specialize the values of the remaining free spins. 

To conclude this Section, we compare briefly our calculation with the 4-dipole in the SU(2) BF case. Instead 
of lOj-symbols, the vertex weight is a 15j-symbol (one additional degree of freedom per tetrahedron). As this is a 
topological case, three delta functions are redundant in the group integral formulation, exactly like in our calculation. 
The amplitude can thus be regularized the same way, by fixing the spins on three internal faces. Once all delta 
functions arc taken into account, the special solutions ([77| holds at the vertices A and B. Therefore, there is no 
integral like in the right hand side of (|M)) . This integral is really the remnant of the way the BF theory is modified to 
get the BC model (i.e. imposing the simplicity constraints in a specific way), which survives even when the amplitude 
is restricted by hand to the special BF-like solutions ((77|) . 



2. Using the tetrahedral graph 



We consider a piece of triangulation formed by four 4-simphces, labeled 1,2,3,4, connected to one another. The 
tetrahedra they share (called internal) are therefore labeled by pairs (jj), 1 < i < j < 4. The boundary has eight 
tetrahedra, each 4-simplex contributing to two, denoted iA^ iB. The two boundary tetrahedra of the simplex i share 
a triangle labeled {AiB), for i = 1,2,3,4. The boundary tetrahedra of the simplices i,j share two triangles, one 
belonging to the tetrahedra of type- A and one to the tetrahedra of type-i?. We label these triangles (ijA) and (ijB) 
and notice that they also belong to the internal tetrahedra (ij). The internal structure has four triangles, which are 
all shared by three 4-simpliccs, and are therefore labeled (ijk), ioi 1 < i < j < k < 4. The triangle {ijk) belongs to 
the three internal tetrahedra (ij), (jk), (ik). 

In the dual picture, 4-simplices are vertices, tetrahedra lines and triangles faces. The 2-complex, denoted r4, is 
depicted (as a graph) in the Figure [5J It has four vertices i = 1,2,3,4, connected to one another by six (internal) 
lines (ij), for 1 < « < j < 4. The tetrahedra on the boundary of the gluing are represented by eight half-lines labeled 
{iA), {iB) for i = 1, 2, 3, 4, each vertex i having two of them. The external faces are broken faces dual to the boundary 
triangles. The external face (AiB) goes along the half-lines (iA) and (iB), and there are four of them. There are 
twelve other external faces, labeled (ijA) (six of them), and (ijB) (six others), for 1 < i < j < 4. The face (ijA) 
goes along the half-line (iA), then the internal line (ij) which connects the vertices i to j, and continues along the 
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Figure 2. There are four 4-siinplices (represented as vertices), all connected to one another by a tetrahedron (represented as 
a edge). 



half-line (jA) (similarly for (ijB)). The internal faces all have three vertices, which allows to have them labeled {ijk) 
and there are four of them. 

To write the spin foam amplitude on r4, we fix the spins of the external faces: ji for the face (AiB), jijA and jijB 
for the faces {ijA), (ijB). The internal faces are closed with three lines, and thus carry effective weights of faces with 
n = 3. The line (ij) has two group elements hij, hji, respectively associated to the half-line connected to i and to j. 
The external half- lines have group elements hi a, hiB- The partition function is 

4 

\[ dh,Adh,B n dh^J Yi [dlk{^J)dJJ(^k^) 5{^^jKk 7fe(zj) hk^hlj 7-i^.j 7j(H) hji^hjl 7"!^^) 

4=1 i^j l<i<j<k<i 

4 

WxiAKl^iB) n Xj,jAihiAh~j^)xj,jA{lT'j^h~l) X3^,BihiBh~j^)xj^,Bih]ihjB). (97) 

i=l l<i<i<4 

Notice that the internal faces form a spherical bubble, identical to the boundary of a tetrahedron. Therefore, the 
contribution of this bubble can be evaluated as an application of the result of the Section IIVI We ignore the external 
faces (putting their spins to zero), and use the formula (|40|) for the BC model on a surface of Euler characteristic 
X = 2 with E = 6 hnes, to get 



Zr,(0,0,0) = C(8) = 



9450 



(98) 



This is obviously finite, meaning that the four deltas in (|97p are all independent. However, if we can project onto the 
special configurations (fTT]). there would be a gauge symmetry of the BF type, which in the case of r4 corresponds to 
one redundant delta (like for any spherical bubble in the BF model). To project onto solutions of the form ([77| . we 
proceed like in the 4-dipole case. We change some of the characters of the external faces with deltas. 
Let us consider 



f ^ ^ - 

•^ 1=1 i^j l<i<j<k<4 i=l l<i<j<4 



XniA ih3AhJ) XuiA (hiAh^l) Xj42a (hiAhJ) 



n n ^i^^^hr^'^ 



i=l,2j"=l,2,3,4 



Sih^Ah^i) 5{hzAh^l) 5{hiAh^l) 



5\hi2h-^^ 73(12)^31^32^3(12) T2(31) /«23 ^2/72(31)) ^y^lihlA 74(13)^41/^43^4(13) 73(41)^34/131^3(41) 
5\hi2h^l 74(12) hiih'^^ 747(12) 72(41) /l24/l21^ 72(41))Xj234 [h23h^l 74(23) /i42/l43^ 74(23) ^3(42) h^^h^^ 73^(421 ) ' i^^) 



(42) 

The group variables we integrate are the same as in (P7| and the products of these group elements appearing in the 
integrand are also the same. Only the functions differ. The four external faces (AiB) are untouched, as well as the six 
external faces {ijB). The characters Xji-A along the wedges of the external faces {ijA) have almost all been replaced 
with deltas, except for the wedge of the face {23A) at the vertex 3, the wedge of the face (14^) at the vertex 4, and 
the wedge of the face {2AA) at the vertex 4. Finally, the delta on the internal face (234) has been changed with a 
single mode Xj23.n to avoid a divergence due to a gauge symmetry as we will see. 
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Now let us consider the effects of all these new deltas in Ip^- We have hi2 = /113, ft.21 
constraint on the face (123) simplifies to /131 = /132, 



/i23- Therefore, the 



S{hiAhj^2) ^i^T-lAhJ) S{h2Ah2i)S{h2Ah23) Si^hi2hj^i 73(12)^31^32^3(12) 72(31) /i23/i2lS2(31) 

== SihiAK^') 5{h^AK^) 5{h2Ah2i) 5{h2Ah2^) 5{h:iih^^). (100) 



Similarly, we have hi^ = /ii4, h^i = /i34 which means that the constraint on the face (134) simplifies to h^i = ^143, 



5{hiAh-i^^) 5{hiAhil) 6{h^Ah^l) SQl^Ah^l) Sl^huh^l 74(13)^41^-43^4(13) 73(41) /i34/i3iS3(4i)^ 

= 6{hiAKi) 6{hiAKl) 6{hzAh^l) SihsAh^l) 5{hiihll), (101) 



and for the face (124), we get /i42 = /141, 

5{hiAh'{^) 5{hiAh'{l) (5(/l2A^^/) '5(/i2A^M^) ^y^l2hil 74(12)^41^12^^(12) 72(41) /i24^^iS^(4i)) 

= 5{h^AK^) 6{hiAKl) 5{h2Ah2l) 5{h2Ah2l) 5{hi^hl^). (102) 

As a result of all the deltas in ([M|) . we find the set of solutions of the contraints, 



^Ti = {hi2 = h 



13 — "■14 



hlAi h2l — h23 — /i24 — h2A, ^31 — ^32 — /l34 = ^3A, ^41 = /i42 = /l43 = /i4A} ■ (103) 



and the 7s can take arbitrary values. The character on the fourth face, (234), is thus simply evaluated on the identity. 



Xj234 \h2-ih^l 74(23) ^42/^43^ 74(23) 73(42) h^^^h^^ ^3(i2) 



Fr4 



(104) 



Clearly, a delta on that face would have been redundant, and caused a divergence of the type 5(1) = ^ ^?234' This 
is the sign of the gauge symmetry which exists when projecting onto Ty^ ■ 

The other characters going along the half-lines (iA) simplify. 



Xi32A(''-3A^32 )Xi41A(^4A^41 )XJ42a(^4A^42 )|JE'r4 — '^J32A '^JAIA dj42A- 



(105) 



The only remaining non-trivial part is the integrals over hiB- Performing the changes of variables hi = hish^ , it 
finally becomes 



^r4 "J234 ^J32A ^J41A "'ji2A 



n 



3¥=i 



(106) 



where we recognize these integrals as the same as the ones on the right hand side of 

A second way to evaluate /r4 is by expanding all deltas as (5 = ^ djXj and integrating the group elements 7 using 
the orthogonality relation ()A4p . For each vertex i = 1,2, 3, 4, we get a lOj-symbol, 



/ dh^Adh,B Yldhij Xj,{h,AhiB)Y[xj,,A{h^Ah^/)x3,,B{hiBhi/) Y[ Xj^.Jhijh^k 



(107) 



j<k 
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Equating the two ways to evaluate Ir^ finally gives 



^J234 "'J32A "-jilA "'ji2A J_ J_ 



/ dh, Xnif^i) Y[xn,B{ht) 



(P. 



'r- E 



J234 



J123 J134 J124 



/ ^ ^J12A^J13a'^J14A 



/ J ^J21A"'J23a'^, 



]12A,J13A 
jliA 



/ , "'331 A ", 



J31-4J34A 




il23 '^il24 '^jl34 



'33AA 



'34.3A 




(108) 



In summary, as in the case of classical discrete gravity, there are special solutions to the delta function constraints 
describing the BC model around which gauge symmetries occur. By changing the partition functions appropriately 
we can enforce a projection onto these special solutions. These altered partition functions will have gauge symmetries, 
and this can be used to derive equations involving the vertex amplitude of the BC model. 



VIII. DISCUSSION 



We discussed possible divergences in spin foam models, in particular in the Barrett-Crane model. The reason for 
considering this model is that the relation of divergences to gauge symmetries is easily found in such a model, for 
which the partition function can be rewritten as an integral over a space of flat connections. This also explains our 
choice of edge weight factors and face weights. 

We presented a simple method to determine the (single bubble) divergences for general spin foam models. We 
argue that to this end one just needs to evaluate the spin foam model on a two-dimensional surface and we gave an 
explicit formula for a large class of models, encompassing BC and EPRL/FK. This allows to deduce the influence of 
the edge and face weights on the convergence of the models very easily. We noticed that requiring invariance under 
face divisions and in addition under (two-valent) edge subdivisions leads to triangulation invariant 2D models, which 
have however divergent spherical bubbles. 

We discussed in detail the effective face weights as they capture the basic, possibly distributional, building blocks 
for the models. We found that for the BC model the effective face weights (for faces with more than two edges) are 
(almost everywhere) regular functions. For the EPRL/FK models, finiteness depends on the choice of edge and face 
weight factors. However applying our arguments on how to evaluate spherical bubbles, one can consider the square 
of the effective face weights. This will identify a distributional character of the effective face weights also for these 
models. The (Lorentzian) EPRL/FK model has been argued to be finite for the 1-5 move and have only a logarithmic 
divergence for the 4-dipole [10| . This suggest that the effective face weights for faces with more than two edges might 
also be finite functions and not be distributional in this case. 

For spin foams recent work [33[ suggested the notion of wave front sets, which specifies the non-smooth part 
of a distribution, in order to study the large spin limit and regularization issues. The wave front sets have been 
identified for both the Barrett-Crane and the EPRL model in [SJl, thus there is clearly the potential for divergent 
behavior. Indeed if the wave front sets would correspond to divergences, regularization is needed (to define products 
of distributions). In that case, the conclusions of [33| for the large spin limit regime would have to be reconsidered, 
as this work assumed that products can be formed from the distributions that occur as amplitudes in spin foams. 
Here we found that (for higher than two-valent faces) the wave front sets for BC correspond to non-continuous or 
non-smooth but also non-divergent behavior (modulo sets of measure zero). 
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Furthermore we analyzed in detail the 4-dipolc case for the BC model. With our choice of edge and face weight 
the spin picture does not allow for a definite conclusion regarding finiteness. We therefore used the fact that the 
BC amplitudes can be rewritten as integrals over some sets of flat connections, which renders the problem accessible 
through the method of [131 ■ We showed that no redundancies of delta functions arise, thus excluding divergences due 
to this reason (another possible source are the measure zero singularities though). 

Redundancies of delta functions would be the sign of gauge symmetries. As such redundancies are not occurring 
- and in addition we found convincing arguments that the BC model is finite on a regular (i.e. involving only faces 
with more than two faces) triangulation - we have to conclude that gauge symmetries, which could be connected to 
diffeomorphisms, are not present. This even seems to hold for configurations, for instance the 5-1 Paehner move, for 
which the symmetries exist on the classical (Regge) level |19j . 

There are however special solutions (of measure zero) for which delta function redundancies can be identified. 
Similarly there are measure zero solutions in gravity (the flat solutions) around which (linearized) gauge symmetries 
can be found. We discussed those special solutions and described the related gauge symmetries for the BC model. 
We developed a method to derive associated Ward-identity-like equations on the vertex amplitude. This is the first 
proposal which enables to extract constraints from a quantum theory with broken gauge symmetries, i.e. equations 
that have to hold for the boundary wave function (which modulo measure factors can be identified with the vertex 
amplitude), extending this way the tools introduced [22| for topological theories. 

The question arises whether those special symmetries can be the seed for the occurrence of more general symmetries, 
that might emerge under coarse graining j43| . The heuristic argument is that coarse graining leads to an effective 
description of the coarse model on a much finer triangulation. On this fine triangulation the curvature per building 
block is very small, so that one is near the flat case, that is on the special solutions around which gauge symmetries 
do exist. 



This mechanism actually works for classical systems 14J] as well as for ID quantum systems [45[. In this ease the 



amplitudes might become more and more divergent under coarse graining, as is indeed the case for ID diseretized 
quantum systems [43. We will leave this question for future work. 

Another interesting question is whether the wave front analysis performed in |33| can be used to analyze the gauge 
symmetries around special solutions also for more general spin foam models. Wave fronts are a refinement of the 
singular support to the co-tangent space. We also used the co-tangent space in ([7T|) to define the notion of gauge 
symmetries we applied in this work. This might make the methods presented here applicable to other spin foam 
models as well. 

Appendix A: SU(2) calculus 

We parametrize group elements as ft. = e*^ "'^ = cos 6I + i sin Oh-a, where n £ 5*^ is the rotation axis and S [0, tt] 
is the class angle. The vector a = {(Jx,(Ty,(Tz) is the 3- vector formed by the Pauli matricef°l. which transforms as a 
CO- vector under the adjoint action, 

j=x,v,z 

R{g) being the rotation matrix in the vector representation, g S SU(2). From this, the orbit of the adjoint action on 
the group is found, 

ghg-^^ COS6II+ t sine {R{g)fi) ■ a ^ e''^ iRia)n)-3 ^ ^^3^ 

meaning that g rotates the rotation axis of h without changing its class angle. 

The matrix elements in the irreducible representation of spin j G N/2 satisfy the orthogonality relation 



dh D<;i\i^ (h) D^£}n, (h) - ■]- Sj,j, S^m.rn, Sn,n, ■ (A4) 

'2) «ii 

Here dh is the normalized Haar measure, dj = 2j + 1 is the dimension of the representation and D^^^ the Wigner 
matrices. A useful property is _D,„„(ft,^^) = Dnm{h) for any h G SU(2). The character in the representation of spin j 



They read 



? !)' ''y=U 0)' "- = (! -ij- (^1) 
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is the trace, 

It satisfies Xj (I) = ^^i ■ The convolution of characters is 

/ dg XjiiK^g) Xj2(.9"^^2) = Sj,j^ — Xh{hi^h2). (A6) 

Functions in L^(SU(2), dh) admit expansions over the Wigner matrices, f{h) = Yl,j£m/2 Em n=-i \/dj fmnDmn{h), 
which is the Fourier expansion. Class functions arc the functions invariant under the adjoint action, so that they only 
depend on the conjugacy class, i.e. the class angle. Characters provide a basis of class functions. The Dirac delta 
over SU(2) is the distribution such that /guc2) dg5{g) f{g) = /(I) and it has the expansion 



%)- E '^^^^(9)- (A7) 

jeN/2 



The delta over the conjugacy class of angle ip is 



S^ih)= / djS{h-fg^-f ^) = E Xji9i>) XjW, (A8) 

where g^, S SU(2) is any representative of the conjugacy class. 
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